Anisotropic induced gravity and inflationary universe 



W. F. Kao 

Institute of Physics, Chiao Tung University, Hsinchu, Taiwan 

Existence and stability analysis of the Kantowski-Sachs type universe in a higher derivative in- 
duced gravity theory is studied in details. Existence of one stable mode and one unstable mode 
is shown to be in favor of the inflationary universe. As a result, the de Sitter background can be 
made to be stable against anisotropic perturbations with proper constraints imposed on the coupling 
constants of the induced gravity model. 
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I. INTRODUCTION 

\o ■ 
o ■ 

Inflationary theory provides a natural resolution for the the flatness, monopole, and horizon problems of our present 
£NJ , universe described by the standard big bang cosmology [1]. In particular, our universe is homogeneous and isotropic 
£j ' to a very high degree ofprecision 0, [3] . Such an universe can be described by the well known Friedmann- Robertson- 
Walker (FRW) metric 0]. 

Moreover, gravitational physics could be different from the standard Einstein models near the Planck scale 
For example, quantum gravity or string corrections could lead to some interesting cosmological applicastions[5j. In 
particular, investigations have been conducted on the possibility of deriving inflation from higher order gravitational 
' corrections @, H, 0, G3- 

, For example, a general analysis on the stability conditions of gravity theories could be useful in screening physical 
models compatible with our physical universe. In particular, the stability condition for a variety of pure gravity 
theories as a potential candidate of inflationary universe in the flat Friedmann-Robertson- Walker (FRW) space is 
derived in Ref. E3, 

In addition, the highly isotropic universe should evolve from some initially anisotropic state before it becomes 
isotropic to such a high degree of precision. Nonetheless, it is interesting by itself to study the stability analysis 
of the anisotropic space during the post-inflationary epoch even anisotropy can be smoothed out by the proposed 
, inflationary process. One would like to know whether our universe can evolve from certain anisotropic universe to a 
stable and isotropic final state. In particular, it is known that such inflationary solution exists for an NS-NS model 
with a metric, a dilaton, and an axion field [111 ]. This inflationary solution is also shown to be stable against small 
perturbations [T3]- Note that similar stability analysis has also been studied for a various of interesting models [H, [r|. 

The importance of higher derivative models have been known in many aspects. In particular, higher derivative terms 
are known to be important for the Planck scale physics fiol. [l5j. For example,_higher order corrections from quantum 
gravity or string theory have been considered as possible inflationary models [171 ] . In addition, higher derivative terms 
also arise as the quantum corrections of the matter fields [13] • Therefore, it is important to study the implications of 
5_| ' the stability analysis of all possible higher derivative models. 

Recently, there are also growing interests in the study of Kantowski-Sachs (KS) type anisotropic spaces (IH fl9l I20I] . 
We will hence try to study the problem of existence and stability conditions of an inflationary de Sitter final state 
for some higher derivative model in Kantowski-Sachs spaces. In particular, a large class of pure gravity models with 
inflationary KS/FRW solutions was presented in Ref. j21j. Any KS type solution leading itself to an asymptotic FRW 
final state will be referred to as the KS /FRW solution in this paper for convenience. 

It has been shown that the existence of a stable de Sitter background is closely related to the choices of the coupling 
constants. Indeed, a pure gravity model given below: [2l| 

C = -R- aR 2 - miR^ + iR^ R P l P R°% (1) 

admits an inflationary solution with a constant Hubble parameter determined by Hq — l/4y if 7 > 0. Here a, (3, and 
7 are coupling constants. This shows that (a) the 7 i?-cubic term determines the scale of the inflation characterized 
by the Hubble parameter Ho and (b) quadratic terms are irrelevant to the scale Hq in the de Sitter phase. The 
quadratic terms are, however, important to the stability of the de Sitter phase. 
Indeed, perturbing the KS type metric with H L — > Hq + 5Hi, one can show that 

SH t = a exp[^^(l + ft)] + di exp[^^(l - 6,)] (2) 
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for 

Si = y / l + 8/[27-9(6a + 2/3)iJ 2 ] (3) 

and some arbitrary constants C;, di to be determined by the initial perturbations. Here Hi = hi/ai with cii(t) the scale 
factor in i-direction. We will describe the notation shortly in section II. It is easy to see that any small perturbation 
SHi will be stable against the de Sitter background if both modes characterized by the exponents 

A± = -[3fr t/2][i±ai] (4) 

are all negative. This will happen if 5\ < 1. In such case, the inflationary de Sitter space will remain a stable 
background as the universe evolves. More specifically, the constraint (3Hq > 35/18 is required for both modes to be 
stable. 

It can be shown that the stability equation for the anisotropic KS space and the stability equation for the isotropic 
FRW space in the presence of the same inflationary de Sitter background turns out to be identical [1(1 [H, 0- 
Therefore, the stability of isotropic perturbations also ensures the stability of the anisotropic perturbations. The 
stability of the isotropic perturbations for the FRW space is important for any physical models. Unfortunately, 
inflationary models that are stable against any isotropic perturbations will have problem with the graceful exit 
process. Therefore, the pure gravity model may have troubles dealing with the stability and exit mechanism all 
together. 

Instead of the pure gravity theory, a slow rollover scalar field may help resolving this problem. An inflationary de 
Sitter solution in a scalar-tensor model is expected to have one stable mode (against the perturbation in SHi direction) 
and one unstable mode (against the perturbation in 5<j) direction). As a result, the inflationary era will come to an 
end once the unstable mode takes over after a brief period of inflationary expansion. Therefore, we propose to study 
the effect of such theory. 

In particular, we will show in this paper that the roles played by the higher derivative terms are dramatically 
different in the inflationary phase of our physical universe in both pure gravity theory and scalar-tensor theory. First 
of all, third order term will be shown to determine the expansion rate H for the inflationary de Sitter space. The 
quadratic terms will be shown to have nothing to do with the expansion rate of the background de Sitter space. They 
will however affect the stability condition of the de Sitter phase. Their roles played in the existence and stability 
condition of the evolution of the de Sitter space are dramatically different. 



II. NON-REDUNDANT FIELD EQUATION AND BIANCHI IDENTITY IN KS SPACE 



Given the metric of the following form: 

ds 2 = -dt 2 + c 2 (t)dr 2 + a 2 (t)(d 2 6 + f 2 (d)dip 2 ) (5) 

with f(9) = (9 1 sinh 9, sin#) denoting the flat, open and close anisotropic space known as Kantowski-Sachs type 
anisotropic spaces. More specifically, Bianchi I (BI), III(BIII), and Kantowski-Sachs (KS) space corresponds to the 
flat, open and closed model respectively. This metric can be rewritten as 

dr 2 

ds 2 = -dt 2 + a 2 {t){ — T + r 2 d6 2 ) + a 2 {t)dz 2 (6) 

1 — kr z 

with r, 9, and z read as the polar coordinates and z coordinate for convenience and for easier comparison with the 
FRW metric. Note that k = 0, 1, —1 stands for the flat, open and closed universes similar to the FRW space. 

Writing i? M „ = G^ v — T^, Einstein equation can be written as D^H^ — incorporating the Bianchi identity 
D^G^ V — and the energy momentum conservation D^T^ — 0. Here G MI/ and represent the Einstein tensor 
and the energy momentum tensor coupled to the system respectively. With the metric ([6]), it can be shown that the 
r component of the equation D^H^ V — implies that 

H\ = H d e . (7) 

This result also says that any matter coupled to the system has the symmetric property T r r = T e . In addition, the 
equations D^H^ = and D^H^ Z = both vanish identically for all kinds of energy momentum tensors. More 
interesting information comes from the t component of this equation. It says: 



(dt + 3H)H\ = 2H 1 H r r + H Z H Z Z . 
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This equation implies that (i) H\ = implies that H r r = H z z = and (ii) H r r = H z z = only implies (d t +3H)H\ — 
instead of H\ — 0. Case (ii) can be solved to give H\ = constant x exp[— a 2 a z ] which approaches zero when a 2 a z — > oo. 
For the anisotropic KS spaces, the metric contains two independent variables a and a z . The Einstein field equations 
have, however, three non-vanishing components: H\ — 0, H r r = H e 8 = and H z z = 0. The Bianchi identity implies 
that the tt component is not redundant and will hence be retained for complete analysis. Ignoring either one of the 
rr or zz components will not affect the final result of the system. In short, the Hi — equation, known as the 
generalized Friedman equation, is a non-redundant field equation as compared to the H£ = and H z = equations. 

In addition, restoring the gu component b 2 (t) = 1/B\ will be helpful in deriving the non-redundant field equation 
associated with Gu that will be shown shortly. More specifically, the generalized KS metric will be written as: 

dv 2 

ds 2 = -b 2 {t)dt 2 + a 2 (t)(—— + r 2 d0 2 ) + a 2 z (t)dz 2 . (9) 
1 — kr A 

In principle, the Lagrangian of the system can be reduced from a functional of the metric g^ Vl C^g^), to a simpler 
function of a(t) and a z {t), namely Lit) = a 2 a z £(5 /J!/ (a(t), a 2 (t))). The equation of motion should be reconstructed 
from the variation of the reduced Lagrangian L(t) with respect to the variable a and a z . The result is, however, 
incomplete because, the variation of a and a z are related to the variation of g rr and g zz respectively. The field 
equation from varying g tt can not be derived without restoring the variable b(t) in advance. This is the motivation 
to introduce the metric ([9]) such that the reduced Lagrangian L(t) = 6a 2 a 2 £(t/ M1/ (6(t), a(t), a 2 (t))) retains the non- 
redundant information of the H\ = equation. Non-redundant Friedmann equation can be reproduced resetting 6=1 
after the variation of b(t) has been done. 

After some algebra, all non-vanishing components of the curvature tensor can be computed: [2l| 



R u tJ = [-B 1 H i + B 1 {H i +H^)]8), (10) 
R l \ l - BxHiHj e^ m e klm + ^ Jz e k iz (11) 



with Hi = (a/a,a/a,d z /a z ) = (H±,H2 = H\,H Z ) for r, 9, and z component respectively. 
Given a Lagrangian L = ^/gC = L(b(t), a((t), a z (t)), it can be shown that 



L = ^L{R ti tp R i \ l ) = ^L{H h H i ^) (12) 



The variational equations for this action can be shown to be: [2l| 

C + H t {^- t +3H)L l = HtLt + HiL* (13) 
r + (|+3H)»L* = (|+3F)L. (14) 

Here Li = SC/SHi, L l = 8C/8Hi, and 3H = Hi. For simplicity, we will write C as L from now on in this paper. 
As a result, the field equations can be written in a more comprehensive form: 

DL = L + Htt— + 3H)L i - HiLi- HiL 1 = (15) 
dt 

D Z L = L + (— + 3H) 2 L Z - (— + 3H)L Z = (16) 
dt dt 



III. HIGHER DERIVATIVE INDUCED GRAVITY MODEL 



In this section, we will study the higher derivative induced gravity model: 

L = -\<?R - aR 2 - 13R£R; + Jf#% 7 R P l P ~ \ d ^4> ~ V(cj>) = ^ 2 L + L 2 + ^L 3 + L (17) 

with L = -R, L 2 = -aR 2 - flR^R 1 ^, L 3 = R*" R?\ p R^ v and L = -\d^d^<j) - V(<f>) denoting the lowest 
order curvature coupling, the higher order terms, and the scalar field Lagrangian respectively. Induced gravity models 
assume that all dimensionful parameters or coupling constants are induced by a proper choice of dynamical fields. For 
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example, the gravitational constant is replaced by 8irG = 2/(e(j> 2 ) as a dynamical field. In addition, the cosmological 
constant becomes V(<p) in this model. There is no need for any induced parameters for the quadratic terms R 2 
and R^ v because the coupling constants a and (3 are both dimensionless by itself. This action of the system is also 
invariant under the global scale transformation: g^„ — * hr 2 g^ v and cj> — > A(f> with arbitrary constant parameter A. 
The corresponding Lagrangian can be shown to be: 

L = e(/) 2 (2A + B + 2C + D)-Aa [AA 2 + B 2 + AC 2 + D 2 + AAB + 8AC + AAD + ABC + 2BD + ACD] 
-2(5 [3A 2 + B 2 + 3C 2 + D 2 + 2AB + 2AC + 2AD + 2BC + 2CD] + +8-^ [2A 3 + B 3 + 2C 3 + D 3 ] 

+l^-V(4>) (18) 

with 

A = H 1 + Hi, (19) 

B = H 2 + ^, (20) 

or 

C = H 1 H Z , (21) 

D = H Z + H 2 . (22) 

This Lagrangian can be shown to reproduce the de Sitter models when we set Hi — > Hq in the isotropic limit. The 
Friedmann equation reads: 



\t4> 2 DL G + DL 2 + ^DL 3 + e<M>HiL' - 2^^i| - \j> 2 + V{<j>) (23) 
2 tp z <p a 2 

for the induced gravity model. In addition, the scalar field equation can be shown to be: 



4 + 3H 4> + V = e<£L -2^L 3 . (24) 

The leading order de Sitter solution with <j) — 4>o an d Hi = Hq for all directions can be shown to be: 

V = V(<j> ) = 3e'<j> 2 Hl (25) 
V'{<h) = 12e>o^o (26) 

under the slow roll-over assumption <f> « V and Hq4> << V. Here e' = e[l — 8jHq / (c^q)]. Indeed, It can be shown 
that in the de Sitter inflationary phase, the ignored part of the scalar field equation evolves as <j> + 3Hq4> ~ 0. This 
equation leads to the approximated solution 

<P ~ <Po + 4k I 1 - ex P (-3 J ff i)] (27) 

during the de Sitter phase Hi = Hq. This result is clearly consistent with the slow roll-over assumption we just made. 
In summary, the leading order equations give us a few constraints on the field parameters: 

dV 

An appropriate effective spontaneously symmetry breaking potential V of the following form 

V(4>) = ^{4? - M 2 + Ge'H 2 ^ 2 - + 3e'H$cj> 2 (29) 



with arbitrary coupling constant A can be shown to be a good candidate satisfying all the scaling conditions (I28|) . 

The value of Hq can be chosen to induce enough inflation for a brief moment as long as the slow rollover scalar 
field remains close to the initial state <f> = 4>q. The de Sitter phase will hence remain valid and drive the inflationary 
process for a brief moment determined by the decaying speed of the scalar field. The stability conditions associated 
with this effective potential implied by the field equations will be studied in the following section. 
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Note that the local extremum of this effective potential can be shown to be (j> = (local maximum) and <fi 2 — <f) 2 n — 
<Pq — I2e' Hq I 'A < (f>Q (local minimum). In addition the minimum value of the effective potential can be shown to be 

V m =V - 36e' 2 ff 4 /A < V . (30) 

The constraint V m > implies that A0q > 12e'H^. Or equivalently, it implies that § 2 m > 0. In addition, we will set 
e<j> 2 n /2 = l/(8irG) — 1 in Planck unit later in this paper. 

When the scalar field settles down to the local minimum cj) m of the effective potential at large time in the post 
inflationary era, it will oscillate around the local minimum and kick off the reheating process. The scalar field will 
eventually become a constant background field with a small cosmological constant V m = V{<p m )- 

IV. STABILITY OF HIGHER DERIVATIVE INFLATIONARY SOLUTION 

Our universe could start out anisotropic and evolves to the present highly isotropic state in the post inflationary 
era. Therefore, a stable KS/FRW solution is necessary for any physical model of our universe. 

Given an effective action of the sort described by Eq. (f29|) , the scaling constraints (|28j) is required for the existence of 
the de Sitter solution Hi = Hq and the static condition <f> — 4>q. In addition to these constraints, small perturbations, 
Hi = Hio + 6 Hi and <fi — 4>a + <5</>, against the background de Sitter solution (i/jOi^o) ma y a ls° put a few more 
constraint to the stability requirement of the system. This perturbation will enable one to understand whether the 
background solution is stable or not. In particular, it is interesting to learn whether a KS — > FRW (KS/FRW) type 
evolutionary solution is stable or not. 

It can be shown that the perturbation equation for the Bianchi models are identical to the perturbation equation 
for the FRW models. Therefore, any inflationary solutions with a stable mode and an unstable mode will provide a 
natural way for the inflationary universe to exit the inflationary phase. Such models will, however, also be unstable 
against the anisotropic perturbations. Therefore, any inflationary solutions with a stable mode and an unstable mode 
is also a negative result to our search for a stable and isotropic inflationary model. As a result, such solution will be 
harmful for an anisotropic space to reach an isotropic FRW space once the de Sitter phase start to collapse. It will 
be shown shortly that the higher derivative induced gravity theory could hopefully resolve this problem all together. 

In practice, perturbing the background de Sitter solution along the SHi direction should be stable for at least a brief 
moment of the order AT ~ 60iT ~ 1 for a physical model. As a result around 60 e-fold inflation can be induced before 
the de Sitter phase collapses. And the resulting universe is stable against isotropic and anisotropic perturbations. In 
addition, the scalar field is expected to roll slowly from the initial state <j> — (f>o- Therefore, the perturbation along the 
6(j> direction is expected to be unstable. This will favor the system for a natural mechanism for graceful exit. Hence, 
we will try to study the stability equations of the system for small perturbations against the de Sitter background 
solutions. 

The first order perturbation equation for DL, with Hi — > Hq + SHi, can be shown to be: 

5{DL) = < HiL ij 5Hj > +3H < > +5 < H,V > +3H < (H.L) + L J )5H,j > 

+ < HiV > 5(3H)- < HiLijSHj > (31) 

for any DL defined by Eq. (fT5]) with all functions of Hi evaluated at some FRW background with Hi = H - The 
notation < AiBi >= Jw=i z AiBi denotes the summation over i = 1 and z for repeated indices. Note that we have 
absorbed the information of i = 2 into i — 1 since they contributes equally to the field equations in the KS type 
spaces. In addition, Lj = S 2 L/SHiSHj and similarly for and L 1 ^ with upper index 1 and lower index j denoting 

variation with respect to Hi and Hj respectively for convenience. In addition, perturbing Eq. (|16p can also be shown 
to reproduce the Eq. ([31]) in the de Sitter phase [l|. 
In addition, it can be shown that 



< HiL a > = 2 < HiL iz >, (32) 

< HiL\ > = 2 < HiL\ >, (33) 

L 1 = 2L\ (34) 

< H.La > = 2 < H,L lz >, (35) 



in the inflationary de Sitter background with Hq — constant. Therefore, the stability equations (|3ip can be greatly 
simplified. For convenience, we will define the operator T>l as 



V L 6H =< HiL a > 8H + 3H < H l L il > SH + 3H < HiL\ + L 1 > 5H + 2 < H t V > 5H- < HiL iX > SH = 0. (36) 
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As a result, the stability equation (|3Tj) becomes 



8{DL) = V L {5Hi + SH z /2) = -V L {SH) = 



(37) 



with H — (2Hi + H z )/3 as the average of Hi. 

Hence the leading order perturbation equation in 8H and 5<f) for the Friedmann equation of this model can be 
shown to be: 



7#o 4 



'yH 4 

6e(l-24l-f)0oHo[^-^o' s 



^S(DL )+S(DL 2 ) + ^S(DL 3 ) = ^ e <j> 2 V SH+^D 2 SH+^V 3 SH (38) 



with VqSH = T>l 5H, V28H = T>l 2 8H and U3SH = T>l 3 8H as short-handed notations. This equation can further 
be shown to be: 



e(l - 24^-£)<&)N> ~ H 8(j)] = 4(3a + /3 - 6 7 -§-)(<5if + 3ff <Sff) + (24 7 -§- - ec^)<5i? 



(39) 



Similarly, the leading perturbation of the scalar field equation can be shown to be: 



3H 5<j) + (V" - \2eHl - 384 7 -£-)<5</> = 6e'(p (8H + AH 5H) 



(40) 



The variational equation of a z can be shown explicitly to be redundant in the limit Hi = H + 5Hi and <j) = 4>q + 5(f) 
following the Bianchi identity. 

Assuming that 6H — exp[hHot]5Ho and Scj) = exp[pHot]5(f)o for some constants h and p, one can write above 
equations as: 



e(l-24^ [p-W = 4(3a +ft - 6^)H [h 2 + 3h + J* 1 ** 1 ** J$#M (41) 
e(j>l % 4(3a + f3 - 67^/^) 



H 



P 2 + 3 P+S -12e'-384 7 |f 



6e'(f> [h + 4:]8H. 



(42) 



These equations are consistent when all coefficients vanish simultaneously. This implies that h = —4 and p = 1. This 
set of solution (h,p) = (—4, 1) hence imposes two additional constraint 



e - 16(3a + /3)§ + 72 7 ^ = 0, 



(43) 



A = " ^ (44) 

with 2A0g = F " - 12e'ij2. 

The coupling constant A has to positive in order for the effective potential V(<p) to be free from run-away negative 
global minimum at cf> — > 00. As a result, the constraints e > and A > imply that 

Therefore, the inflationary phase will remain stable against small perturbation along the SH(= exp[— AHqO^SHq) 
direction. In addition, the inflationary phase also has an unstable mode when we perturb the system along the 
5(j){= exp[.ffot]<50o) direction. The unstable mode is in fact consistent with the slow rollover assumption. The scalar 
field is expected to roll slowly off the initial state (f>o for a brief moment during the inflationary era. This unstable 
mode is hence responsible for the graceful exit process. 

Hence such system with a stable mode and an unstable mode is a very nice candidate for a inflationary model. The 
stable mode 8H = exp[— 4Hot]SHo implies that the de Sitter background solution will remain stable against small 
isotropic perturbation. It also implies that the system is stable against any anisotropic perturbation along all 8H L 
directions. Therefore, the induced gravity model with a scalar field can indeed resolve the stability problem of the 
pure gravity model. 
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V. CONCLUSION 

The existence of a stable de Sitter background is closely related to the choices of the coupling constants. The 
pure higher derivative gravity model with quadratic and cubic interactions [2lT | admits an inflationary solution with 
a constant Hubble parameter. Proper choices of the coupling constants allow the de Sitter phase to admit one stable 
mode and one unstable mode for the anisotropic perturbation. 

The stable mode favors a strong inflationary period and the unstable mode provides a natural mechanism for the 
graceful exit process. It is also found that the perturbation against the isotropic FRW background space and the 
perturbation against the anisotropic KS type background space obey the same perturbation equations. This is true 
for both pure and induced gravity models. As a result, the unstable mode in pure gravity model also means that the 
isotropic de Sitter background is unstable against anisotropic perturbations. Therefore, small anisotropic could be 
generated during the de Sitter phase for pure gravity model. 

We have shown that, for induced gravity models, stable mode for perturbations along the anisotropic SHi direc- 
tions does exist with proper constraints imposed on the coupling constants. In addition, another unstable mode for 
perturbation against the scalar field background 4>q also exists with proper constraints. Therefore, isotropy of the de 
Sitter background can remain stable during the inflationary process for the induced gravity models. 

Explicit model with a spontaneously symmetry breaking </> 4 potential is presented as an example. Proper constraints 
arc derived for reference. It is shown that the scaling conditions (28) must hold for the initial de Sitter phase solutions 
in order for the system to admit a constant background solution during the inflationary phase. Additional constraints 
(|43|) and (|44|) are also shown explicitly for this model. In addition, the inequality shown in Eq. (|45|) is expected to 
hold for a physical model with this potential. 

In summary, we have shown that a stable mode for (an)isotropic perturbation against the de Sitter background 
does exist for induced gravity model The problem of graceful exit can rely on the unstable mode for the scalar field 
perturbation against the constant phase 4> . 

It is also found that the quadratic terms will not affect the the inflationary solution characterized by the Hubble 
parameter H . These quadratic terms play, however, critical role in the stability of the de Sitter background. In 
addition, it is also interesting to find that their coupling constants a and (3 always show up as a linear combination 
of 3a + (3 in these stability equations. Implications of these constraints deserve more attention for the applications to 
the inflationary models. 
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